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EXPERIMENTS ON THE RECOVERY OF WASTE HEAT 



IN COOLING DUCTS 
By Abe Silverstein 



SUMMARY 



Tests have "been conducted in the N.A.C.A. full-scale 
wind tunnel to investigate the partial recovery of the 
heat energy which is apparently wasted in the cooling of 
aircraft engines. The results indicate that if the radi- 
ator is located in an expanded duct, a part of the energy 
lost in cooling is recovered; however, the energy recovery 
is not of practical importance up to airplane speeds of 
400 miles per hour. Throttling of the duct flow occurs 
with heated radiators and must he considered in designing 
the duct outlets from data obtained with cold radiators 
in the ducts. 



INTRODUCTION 



The possibility of recovering a small part of the 
heat energy wasted in the cooling air of aircraft engines 
has been indicated by theoretical considerations (refer- 
ence .1). Calculations show that if the waste heat is added 
at a pressure above that of the external stream to the air 
flowing through a duct, the momentum of the air is increased 
and a thrust may be derived. Experiments have been made 
in the N.A.C.A. full-scale wind tunnel to verify this ener- 
gy recovery in wing cooling ducts as an extension of the 
previous tests (reference 2) which were made with cold ra- 
diators . 

The wing and duct arrangements are the same as de- 
scribed in reference 2. The heat was added in the duct by 
an electrical resistance unit used to simulate the radia- 
tor. Tests wore made of four of the ducts in which the 
effects of adding heat were determined by measurements of 
the drag, quantity of air flow through the duct, air tem- 
perature rise across the heater, static pressure at the 
heater, and the electrical input. The section drag coef- 
ficients of the duct were measured by means of the momen- 
tum method to avoid the inaccuracy arising from measure- 
ments of small changes in the duct drag as differences in 
the drag on the complete wing. 
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The experiments were of necessity planned to utilize 
the available equipment and power resources, and optimum 
tunnel speeds and power input could not "be attained. The 
tests should therefore "be considered as preliminary to in- 
dicate the nature of the effects and to develop the test- 
ing technique. 



SYMBOLS 

Q, heat content or mass flow. 

7, ratio of specific heats, Cp/C y . 

Cp, specific heat at constant pressure. 

T, absolute temperature. 

p, absolute pressure. 

t|q, efficiency of heat cycle. 

T|jj, efficiency of momentum recovery. 

Tij , over-all efficiency of heat exchange. 

p, density. 

a, density ratio. 

V, air speed. 

Vg , velocity at the radiator based on the frontal 
area of the radiator (Vg = for cold 

radiator) . 

V x /V, flow ratio. 

q, dynamic pressure. 

c, chord of wing at . cooling duct. 

H, totsl head. 

c 4 , section profile drag coefficient, 
o 

Ag, radiator frontal area. 
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K, ratio of. pressure drop across radiator to dynamic 
pressure at face of radiator (AP/qjj). 

P, power input. 

AP, pressure drop across radiator. 

AF , power recovered from decrease in drag. 

AP power to equalize flow. 

a 

AP , power loss due to heating of radiator. 

3 

AP , power regained from heat cycle. 

RESUME OP THEORY 
Efficiency of Heat Recovery 



If heat is added to the air flowing through a duct 
(fig. 1) the increase in kinetic energy of the air equals 
the heat added in the duct, Q, minus the heat retained 
in the exit, The efficiency of the heat cycle, r\^, 

is therefore equal to •■- ■ ' • 



q 



Since 



. Q = C p (T a - T x ) 



and 

Qi = C p (T 3 - T 0 ) 

then 

/HQ = 



(T a - T x ) - (T 3 - T 0 > 



for the adia*batic change in the expanding and contracting 
sections of the duct, 
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therefore 
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Although this efficiency, r^, indicates the propor- 
tion of the heat energy available to increase the kinetic 
energy of the flow, a further efficiency of conversion of 
the kinetic energy into thrust must be considered. The 
increase in kinetic energy per second is equal to 

l/2 m (V 4 2 - V 2 ). The force on the duct is the increase 
in momentum per second equal to m (V 4 - V ) . Since the 
airplane velocity is V Q , the energy recovery from this 
thrust equals m V 0 (V 4 - V 0 ) . The efficiency of the con- 
version of the kinetic energy into thrust power equals 

- ^ * % = 2 = ^ ' 

^ (; ;;,;v ; ,- ; ^ W 1/2 m (V 4 2 - V Q S ) v 4 + V o 

The over-all efficiency of the heat exchange, T\qT\ m , 
equals 



(2) 



In the usual case, V Q = V , and % is about 1. A 
curve of T\q against V is shown in figure 2. In comput. 
ing the value of T|q for this curve it has been assumed 
that 0.8 of the free stream dynamic pressure is converted 
into static pressure in the duct. 



Power Loss in Heated Radiator 



The net gain realized by adding heat in a duct i..s 
decreased by the higher power consumption in the heated 
radiator. The power absorbed by a radiator is propor- 
tional to PrVr 3 - For equal cooling with a constant tem- 
perature difference at the radiator the PrVr is essen- 
tially independent of the absolute air temperature, and 
since 

T i 

D = P — — 

M R H l 



then 



and 



1 1 



1 



If the mean radiator temperature, Tp, , is assumed to be 
T 1 + T 3 



then 



P V 3 
R R 



T + T x 2 

(3) 



"l 



For a rise in air temperature T g - of 100° F., the 

radiator loss may be increased as much as 20 percent. 
This power loss defined as AP 3 is independent of air 
speed and is, therefore, of greatest importance at the 
lower air speeds in which case the energy recovery from 
the heat cycle is only a fraction of the energy lost in 

/ T + T \ 8 

the radiator. The value of the function f _i 2 j — x 

>> 2T X ' 

which is the ratio of the increased power absorbed in the' 
hot radiator to the power absorbed in the cold radiator 
is shown in figure 3 for various values of T 2 - T 1 . 



Throttling of Flow :in Heated D\icts 

The addition of heat in a radiator duct also serves 
to throttle the air flow. It may he readily shown that 
adding heat to' a slow-moving air stream does not greatly 
change the total pressure, so that at the outlet the total 
pressure is ahout the same with either the hot of the cold 
radiator. The static pressure at the outlet is independ- 
ent of the conditions within the duct, so that the dynamic 
pressure at the outlet must also "be the same for either 
case. (Subscripts h and c in the following equations 
refer to hot and cold radiators.) 

Since 
and 

T 3 

n 3 c J 

3 h • 

then 



7 ^ 



therefore 



P 3 V 3 

3 h 3 h 



and 



Ap 3 v 3= i _ 



For -a temperature rise of the heated air of ahout 
100° F. the mass flow will he decreased ahout 10 percent. 
The outlet area of wing cooling ducts should therefore he 
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slightly larger than indicated "by tests with cold radia- 
tors. Values of the function 1 - / , which repre- 




sents the decrease in flow due to the heat, are plotted 
in figure 3 for various values of T 2 - T . 

EXPERIMENTAL METHOD AND TESTS 

The tests were conducted in the N.A.C.A. full-scale 
wind tunnel which is described in reference 3. The duct 
and wing arrangement shown in figure 4 are described in 
reference 2, and a drawing of the wing showing the in- 
stallation of the heating unit is shown in figure 5. 

Four ducts were investigated which were similar in 
all respecte-except for the inlet size and position. In 
the nomenclature of reference 2, these ducts are: 

Number Designation 

1 4.6- 2a - 4 - 70 

2 6.0 - 2a - 4 - 70 

3 4.9 - 0 - 4 - 70 

4 6.0-0.-4-70 

The ducts will he hereafter referred to as No. 1, 2, etc. 

The electrical heater used in the duct consisted of 
grids of flat bars connected so as to absorb the available 
power of 100 kw. The resistance units were connected so 
that more of the power was absorbed at the middle of the 
duct than at the edges. 

The temperatures of the air behind the heater were 
measured by means of four calibrated thermocouples. The 
temperature of the air in the wake behind the wing was 
measured by a vertical rack of thermocouples which enabled 
the distribution of heat across the wake to be determined. 

The quantity of air flowing through the duct and the 
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pressure in the rear of the radiator were determined from 
measurements of the total and static pressure over a sec- 
tion of the duct near the outlet. The pressure drop 
through the heater grid was measured for a range of air 
speeds . 

The section drag coefficients of the ducted wing were 
computed from the Jones equation (reference 4) from meas'- 
urements made "by means of a wake-survey rack. The wake 
rack consists of a comb of total head tubes and a separate 
comb of static tubes each about 2 feet long. A simultane- 
ous photographic record of the pressures in all the tubes 
was obtained by means of a multiple manometer. 

The total and static pressure in the duct, air tem- 
peratures behind the radiator, air temperatures in the. 
wake, electrical input to heater, and wake pressures were 
measured for each of the ducts at a lift coefficient of 
C ^ = 0.24 (corresponding approximately to the high-speed 
condition) at a tunnel speed of 83 miles per hour. The 
wake measurements were made for the lateral stations at 
the duct center and 1 foot' from the duct center. 



RESULTS AND DISCUSSION 



In the usual case the section-drag coefficients of a 
wing are obtained "by means of the Jones equation (refer- 
ence 4) as 




( H - p ) 
< H o - Po> 



1 - 



g - PpJ 



(5) 



When heat is added in the wing the air density in the wake 
is changed, and the equation should be modified to 




a(H - p) 
(H 0 - P 0 ) 



1 - 



(H - Pq) 
(H ft - P n ) 



(6) 



in which a equals P. /P . In order to check the assump- 
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tion used in deriving equation (6), measurements of the 
drag and heat content were made at a number of different 
distances behind the trailing edge, for one of the wing 
duct arrangements. These results are shown in figure 6, 
and it may be noted that the drags measured at different 
distances behind the trailing edge are in reasonable 
agreement. The low value obtained at the 72-inch station 
may be due to the lower accuracy of the momentum method 
at the greater distance behind the trailing edge. 

The heat contents at the different stations are in 
good agreement. The effect of the heat in decreasing the 
drag is clearly shown in figure 6 by the low momentum loss 
in the part of the wake in which the heat was concentrated. 

A summary of the test results is shown in table I. 
The reduction in the section-drag coefficient due to add- 
ing heat in the duct is about 0.0030 for all the ducts. 
This drag increment- varied from about 10- to 15 percent of 
the section-drag coefficient of the wing and duct, and' was 
therefore measured with considerable accuracy. 

Representative curves are shown in figure 7 of the 
total and static pressure near the duct outlet, from which 
data the quantity of air flow through the duct was comput- 
ed. The flow through the duct was throttled from 10 to 
15 percent by the addition of the heat. This is slightly 
more than would be calculated from the data of figure 3; 
however, the effect of the increased resistance of the 
heated radiator on the flow which was omitted from the 
previous discussion may readily account for the difference. 
Data showing the effects of varying the resistance in the 
duct on the air flow are shown in reference 2. 

The horsepower input into the duct per foot was ob- 
tained by averaging the integrated value of the heat con- 
tent in the wake at the two measuring stations (fig, 6). 
The temperature difference across the radiator was deter- 
mined from the thermocouples. 

The power recovered per foot of duct width, , was 

calculated from the decrease in drag and the velocity, as 



( A c d q c) V 
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Since 



and 



then 



V = 83.0 miles per hour 



c = 9. 67 feet 



AP ■■= 37.8 Ac, 
l d 

A comparison of the duct resistance, with and without 
heat, must necessarily be made for equal mass flows. The 
comparison is made in table I for the flow corresponding 
to the heated radiator case, and it is therefore necessary 
to credit the cold duct with the difference in internal 
power for the two flows. This power difference, AP 2 , was 
computed on the assumption that the external losses re- 
mained constant for the range over which the V R /v changes. 

The power loss in the radiator is QAP in which AP 
is Kq^. The value of K was determined from measure- 
ments of the pressure drop across the radiator, the values 
of which are shown in figure 8 for a range of air speeds. 

Since 

KPV, 3 At, 
QAP = - 

2 X 550 ' 

then 

AP 3 = KA H P(V ic 3 - 7 Xh> 
1100 

The increased power absorbed by a hot radiator, AP 3 , 
has been discussed in the resume of the theory. Knowing 
the temperature difference, T g - Tj_ , the increased power 
due to heating may be obtained from figure 3 as a fraction 
of the power absorbed by the cold radiator. This incre- 
ment tends to decrease the useful power that may be re- 
covered in the duct. 



In order to compute the power regained due to the 
adiabatic expansion of the air after the heat is added, it 
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is necessary to debit the measured drag differences with 
the power required to adjust the flows to the same quan- 
tity and to credit it with the power lost due to. the high- 
er resistance of the hot radiator. That is, 

AF D = AP, - AP + AP 

The values of Ap^ which are given in table I show 

a considerable scatter which is due largely to the fact 
that they are relatively small differences of larger 
quantities. The cycle efficiency is obtained from 
APq/P = T)q. The experimental r\q obtained from the fore- 
going relation and the theoretical value obtained from 
equation (1) using the pressure behind the radiator for 
the term p x are shown in table I. 

The experimental efficiencies are lower than the the- 
oretical ones and the ratios of the two are somewhat in- 
consistent. As an average of the four cases the experi- 
mental efficiencies are about one-half of the theoretical. 
Lt.will be noted that the efficiencies are low, which is 
due to the relatively slow speed at which the tests were 
made. The efficiency varies as V 3 as shown in figure 2, 

Des ign application. - Prom analysis of the results it 
is indicated that if only the heat from the radiator is 
added in the duct from the standpoint of performance cal- 
culations the effect of the heat may usually be neglected; 
that is, the duct efficiencies obtained with cold radia- 
tors may be applied with sufficient accuracy to represent 
the operating condition. Assuming a high speed of 400 
miles per hour and a conversion of 0.8 of the dynamic 
pressure into static pressure in the duct, a theoretical 
efficiency, T^, of 4 percent is indicated (fig. 2). Pur 
ther assuming the average value of 0.4 for the ratio of 
the heat in the cooling air to the brake horsepower, the 
maximum power return passible is 1.6 percent of the brake 
horsepower. Since the actual efficiency may be only one- 
half of the theoretical, the actual maximum power regain 
may be less than 1 percent of the brake power. From this 
value must be subtracted the power loss due to the heating 
of the radiator and several other minor losses (reference 
1). 

The actual regain in power even for the higher speed 
airplanes may, therefore, be less than 1 percent of the 
engine power which is smaller than the possible accuracy 
of the calculations. 
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The effect of the heat must he considered, however, 
in calculating the quantity of air that will flow through 
the duct. It will "be necessary to provide a larger duct 
outlet and in some cases a larger radiator than indicated 
from tests made with cold radiators. In the present tests 
with temperature differences of ahout 80° to 100° I. 
across the heater, the air-flow quantities were decreased 
from 10 to 15 percent helow those for the cold radiator 
tests . 

This paper has discussed the possibility of recover- 
ing a part of the energy which is apparently lost in the 
cooling of engines. It is of interest to point out that 
if the heat of the exhaust is discharged into the cooling 
duct, a substantial recovery of power may he made. This 
recovery of power will be no greater, however, than that 
resulting from discharging the exhaust gases rearward by 
a suitable exhaust pipe. In effect, the cooling duct be- 
comes an exhaust pipe, and its use for this purpose will 
be primarily dictated by the arrangement and detail design 
of the airplane. For some designs the cooling duct Will 
provide a convenient place for expanding the exhaust gases 
so that they may be discharged more nearly at the free- 
stream velocity. 



Langley Memorial Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. , April 3, 1939. 
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TABLE I 



Duct 


cold 


c d 
hot 


Ac d 


7 1 /? 
cold 


hot 


Pressure 
rise at 
radiator 

Pi " Po 
: lb./ 

sq.ft. 

hot 


Power 
input 
hp./ft. 

P 


Change 
in air 
temp, 
at rad. 

F deg. 


Total 
power 
recov- 
ered 
hp/ft. 
A Pi 


Power 
required 
to pro- 
vide equal 
flow 
hp./ft. 


Power 
lost due 
to heat- 
ing of 
radiator 
hp./ft. 


Power 
regained 
from ex- 
pansion 
hp./ft. 
AP q 


Experi- 
mental 
effi- 
ciency 
of heat 
cycle 


Theoret- 
ical ef- 
ficiency 
of heat 
cycle 


1 


,014? 


.0127 


.0020 


.432 


.398 


11.1 


44,3 


83 


.076 


.094 


,056 


.038 


.086 


.148 


2 


.0154 


.0133 


.0021 


.437 


.383 


10.9 


46.0 


98 


.079 


,127 


.060 


.012 


.026 


.146 


3 


,0180. 


.0159 


.0021 


.354 


.302 


7.7 




54.0 


111 


.079 


.089 


.034 


.024 


.045 


.102 


4 


.0187 


.0164 


.0023 


.331 


.286 


7.2 ! 65.0 


98 


.087 


.069 


.024 


.042 


.065 


.095 
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Air temperature rise in radiator, T 2 -T-j, 
figure 3. 
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Cooling Winc( Plan Form 

showing location and size of heater. 




Section A _ A 

NACA 23017 Centersection shoi/vinq 
installation of heater. 
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Figure 8. 
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